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We complete the derivation of the Cornwall-Jackiw-Tomboulis effective potential for quark propa¬ 
gator at finite temperature and finite quark chemical potential in the real-time formalism of thermal 
held theory and in Landau gauge. In the approximation that the function Aijp^) in inverse quark 
propagator is replaced by unity, by means of the running gauge coupling and the quark mass function 
invariant under the renormalization group in zero temperature Quantum Chromadynamics (QCD), 
we obtain a calculable expression for the thermal effective potential which will be a useful means to 
research chiral phase transition in QCD in the real-time formalism. 
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Effective potential Uiili is a basic means to research the vacuum of quantum held theory and spontaneous 
symmetry breaking. For discussion of chiral wmmetry phase transition in QCD, one must generalize the Cornwall- 
Jackiw-Tomboulis (CJT) effective potential 0 for quark propagator to hnite temperature and finite quark chemical 
potential case. At present, some work on the thermal CJT effective potential are mainly based on the imaginary-time 
formalism 0, very few were found in literature which are based on the real-time formalism |^. In fact, in the real-time 
formalism, besides that it is no longer necessary to do analytic continuation of the discrete Matsubara frequency, one 
can directly use the well-known results in zero temperature QCD, e.g. the running gauge coupling and the running 
quark mass function consistent with the renormalization group (RG) analysis before thermalizing propagators by 
thermal transformation matrices, and this will certainly bring great advantages for practical calculation. 

In a preceding paper for deriving the thermal Schwinger-Dyson (SD) equation for quark self-energy, we gave a 
formal expression for the thermal CJT effective potential for quark propagator in Landau gauge. In this paper, we 
will continue to reduce the formal expression to an explicit formula for future use. 

For convenience of the derivation and self-containment of the paper, we will first put down the relevant results given 
in Ref. 0. The formal expression for the CJT effective potential for quark propagator at finite temperature T and 
finite quark chemical potential fi has been expressed by 

V[G,G*] = Vi[G,G*] + V2[G,G*] (1) 

where in momentum space the contribution from one-loop diagrams 

Vi[G,G*] = ziV/iV,((cos2d4trln[5(p)G-i(p)] +tr[,S-i(p)G(p)]}) 

-b(sin2 6»p{trln [5'*(p)G*"^(p)] -b tr[5'*"^(p)G*(p)] }) - (trl)^, (2) 

and the contribution from two-loop vacuum diagrams 

V2[G,G*] = ^-NfN,G2{NX{9%p-qfMrG}r^{p)YG}r\q)\ [D^.(p - g)] “)) 

-\NfNcG2{Kc){{9‘^[ip - qY]'^T:[l^G}^(p)-i''G‘^^(q)\ [D^,,{p - q)]l^)). (3) 

In the above denotations, (■ • ■) and ((• ■ •)) represent respectively the integration f d^p/(27r)'^ and / d'^pd'^q/{2 tt)^, Nf 
and Nc are respectively the number of flavor and color of the quarks, G 2 {!jLc) = (-^c ~ l)/2fVc, 

Sijp) = i/+ ie) and S*{p) =— ie) with (4) 
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are the propagator of massless quark and its conjugate, 

G{p) = i/\^^ — m{p^) + ie\ and G*(p) = —i/[ ^ — to(p^) — ie] (5) 

are complete quark propagator and its coimigate, in which we have made the approximation AijP') = 1 on the basis 
of the arguments given in Sect.IV of Ref. Q and the replacement B{p'^) —> mljP'), g‘^\{p — g)^] is the running gauge 
coupling, tr means the trace of spinor matrices, 

sin^ 6»p = 6»(/)h(p° - p) + 0{-p°)h{-p° + p) with n{p° - p) = + l] and /? = l/T, (6) 


(p) = cos^ 0pG{p) — sin^ 9pG*{p), G^^ (p) = — cos 0p sin0p + G*(p)] = — e^'^Gy (p), (7) 


and 


[Dp^{p-q)]^ = Df,^{p-q)-gp„2Trn{p°-q°)6[{p-q)'^], [Dp„{p-q)]^ =^ ^/‘^2Trn{p°-q°)6[{p-qf] (8) 
with 

-i . {p- q)iiijp- g)v 


Dpvip-q) = — - ^^2 , . Ap,^{p-q), Ap^^ip - q) = gp^ - ^ 


ip - 0 Y + *£ 


{p - qY + is 


and 


( 9 ) 

( 10 ) 

( 11 ) 

( 12 ) 


Vi^ = *iV/V,(sin2 0p{trln[5*(p)G*-i(p)] + tr[5*-i(p)G*(p)] - tr In [5(p)G-i(p)] - tr[5-i(p)G(p)] }). (13) 

By means of Eqs.(4) and (5) we get 


n(p°-gO) = i/[e/5b“-9”l-l]. 

Now we begin performing detailed calculations of V[G, G*]. First rewrite Vi[G, G*] in Eq.(2) by 

Vi[G,G*] = Vio + ViT 

with 

Vio = iNfN^ (^(|trln [S'(p)G"^(p)] + tr[5'“^(p)G(p)] “ trl)) 

and 


Vio = iNf Nj tr i In 


1 - 


mijP') 


■ le 


+ 


m(jY) 


^ — mijp') + ie 


= iV/V7tri^ 


= i2NfN,J 

= -2NfNc I 


U=i 

(Yp 

w 


-1 

m{pY 

?r 

h 

to 

k 

_p + ie_ 

p — m(jA) + ie 1 


1 In 

nY{pY' 

1 

p^ -\- ie _ 


+ 


2rrY(jP‘^ 


p^ — w?{p^) + ie 


Yp f, p^+TO^(p^) 2m?{p'^) 

In - 


(27r)'^ ( p2+TO^(p2) 

In the last line we have made the Wick rotation and the changes of variables 

/ = ip°, P" =p' {i = 1,2,3). 

In similar way we obtain 


ViT = *2iV/Vc^ sin^ 0p I In 


1 + 


nYijP') 


-2p" 


1 


p2 _ ^2^p2'j _|_ jg. 

1 


-In 


1 


Yp^) 


p2 _ pij2f^p2'j _ jg. 


p 2 _ -\- ie p2 — m?(jp) — ie 

rrY{pY 


p2 _ 


2p^ 

p2 _ pi^2(jj2'^ _|_ 


(14) 

(15) 


= — 4iVyfVc^ sin^ 0plm |ln 

= -ANfNj shY 9p\ - ^nY{pY5[p^ - rrY{pY] V (+ 27rp^S[p^ - m^{pY] 
\ I 'ti\ [p^-m^pYf+e^ ) 

= —47riVyiVc^ sin^ 0prrY{pY5[p^ — m^(p^)] 


( 16 ) 
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where the relation S[p'^ — — TO^(p^)] = 0 has been used. 

Next we rewrite V 2 [G, G*] in Eq.(3) in terms of Eqs.(7)-(9) by 

V 2 [G, G*] = V 20 + V 2 T, 


( 17 ) 


where V 20 and V 2 T are respectively the contributions from the two-loop vacuum diagrams without and with thermal 
corrections of propagators. V 20 can be expressed by 


1 


V 20 = 2( 2) , ■ 2 2(2\,- 

2 \ \ '- J l L J L J J p 2 _ j j _|_ jg 

= -2NfNcC2{Kc){{g^ [(p - <lf] [i'm(p^)m{q^) + {p - qfE{p, g)] 


— w?{p'^) -\- ie q^ — w?{q^) + ie 


(p - qY + *£ 

= -6NfN,C2{NJ 


(Yp(Yq 2r ,_2 -2M rn(p2) 
-g [max(p , g jj ■ 


m{q^) 


(27r)® p2 _|_ 77j2(p2^ g2_|_j^2^g2^ (^p — 

where we have used the result 

trjy'" [ ^ + m(p^)] 7'" [ ^ -f TO(g^)] I - g) = 4[3TO(p^)TO(g2) + {p - qfE(p, g)] 


with the denotation 


E{p, g) = 1 - p^ +q" + 


2 , (p" - <Y? 


(p - qf 


/ip-q)Y 


and after the Wick rotation, the running gauge coupling g^[(p — g)^] has been approximated by 0 

g2[max(p2,g2)] = e{p - q)g'^{p'^) + e{q - p)g^{f) 

SO that 

J dVLqE{p, q) = 0. 

By using Eq.(19) and the result that 

tr| 7 ^[ to(p^)] 7 ‘'[ TO(g^)] = 4 [4m(p^)m(g^) - 2p • g] 


(18) 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


and considering some symmetry of the integrand under the exchange between the variables p and g we may express 
E 2 T by 

V 2 T = i47rlV/lVcG2(iV^)g^(0)^^[4m(p^)m(g^)-p^ - g^]n(p‘’- g°)J[(p - g)^] 

iAttsuY 9p5\p^ — m^(p^)]P 


p^ — m?{p'^) +ie q^ — m?{q^) -\- ie 
9 ^ [(p - qY^ 


g2 _^2(g2) 


-8TTNfN,C2iNY 

I 


— \8m{p^)m{q^) + {p - qYE(j), g)] sin^ 6»p(5[p^ - m^(p^)] 


(p - qY + 

— - TTsin^ 9qS[q'^ - irY{q‘^)] | )) + AE 2 T, 


g2 _ j^2('g2^ _|_ jg 

where ”P” in the second line of Eq.(24) means the principle value of the integration and 

AE 2 T = -il6TT^NfNcC2iNY(^(^g'^[ip-qY] [4m(p^)m(g^) -p^ - q‘^]d[p^ - nrY{p'^)]d[q'^ - rrY{q'^)]d[{p - qf 
n{p^ — g°) |sin^ 9p sin^ 9q — sin^ 9p + sin dp cos dp sind^ cosdqe^l^°“'^''^/^| 


(24) 


( 25 ) 
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Let = m\ is a real root of the equation i.e. mi obeys the equation then we will have 

6[p'^ — m^(p^)] = S{p^ — m\)/f{m\) with /(mf) = |1 — dm'^(jp‘)/dp^\p 2 ^^ 2 . Thus we get 

AV 2 T = -ii2Tr^NfNcC2{Kc)9^{^)j^^{{5{p^ - ml)5{q^ - ml)6[{p - qf] 

n(p^ — g°) |sin^ 0p sin^ 9q — sin^ 9p + sinOp cos 9p sin Oq cos 

= -i32Tr^NfNcC2{NJg^{0)j^^^(^(^S[{p + q)^-ml]S{q^-ml)Sip^) 

n{p^) |sin^ Op+q sin^ 9q — sin^ 9p+q + sin9p+q cos 9p+q sin0q cos 9qe^^^ (26) 

In the second equality above, we have made the change of variables p — q ^ p. It is indicated that the condition in 
which the three <5 functions in Eq.(26) are not equal to zeroes simultaneously is p° = |p] = 0. Taking this and the 
definition of sin^ 9p given by Eq.(6) into account, we are led to that 

[sin^ dp+g sin^ — sin^ 0p+q + sin0p+q cos0p+q sindq cosdqe^l^ IpO-o = sin^ [sin^ 0^ — 1 + cos^ dq] =0. 

On the other hand, it is also noted that because 

J dV(/)^^(p°) = y dp° J dflp-J d\^\p\^^[S{p° - 1^) + d(p° + 1^)] 

no singularity could appear in integrand in Eq.(26) when \p\ 0. Hence we can conclude that 

AE 2 T = 0. (27) 

In summary, from Eqs.(l),(ll),(14),(16),(17),(18) (24) and (27) we may write the final expression for the CJT effective 
potential for quark propagator in the approximation A(jj^) = 1 at finite temperature and finite chemical potential in 
the real-time formalism and in Landau gauge as follows. 

V[G, G*] = Eo + Et (28) 


E) = Eio -I- Eo 

= -2N,nJ 1„ + ) 

\ p^ p^ + m^{p^)/E 

m{p^) 




(29) 


where {■•■)e and {{■■•))e means the integrations of the Euclidean momentums / d‘^p/{2E)‘^ and / d'^pd'^q/{2 tt)^ 
respectively. 

Vt = Eit + Et 

= —AnNfNcl^ sin^ 9pm^{p^)S[p^ — TO^(p^)] ^ 

+i4ENfNcC2iNjg'^{0)(^(^[4:m{p'^)m{q^) - p'^ - g^]n(p° - q°)6[{p - qf] 

- iAesiv? 9p5\p^ — TO^(p^)]P 


p^ — m‘^{p^) +i£ q^ — m?{q^) -|- ie 


q^-m^iq"^) 


-87rA/fVeG2(fy;j<^(| y y y [3TO(p^)w(g^) + {p - qfE{p, g)] sin^ 9p6[p'^ - m‘^{p'^)] 


'’Ai‘ - ’""(■'")]})) 


(30) 


can 


In Eqs.(29) and (30), the running gauge coupling constant [(p —g)^] and the running quark mass function m{p^) vcm 
be taken from the known results in zero temperature QCD which are consistent with RG analysis. After introducing 
the infrared momentum cutoff p^, we can write HD 

^^QCD 


( 31 ) 
















and 
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m{p^) = m In 


\P^\+Pc 

^QCD 


-A /2 


+ 




3 {\p^\ + Pl)e{p^ 


In 


b^l +Pc 

^QCD 


A/ 2-1 


f 2^ _ P 


(32) 


where A = 24/(33 — ^Nf), Aqcd is the RG-invariant mass scale parameter, (j) = (0|7/’'i/'|0)/(ln^^/Ag^^)'^/^ and 
m = m[ln(/i^/AQ(^^)]"^/^ are respectively the RG-invariant quark-antiquark condensates and the current quark mass, 
and /i is the renormalization scale parameter. Eqs. (28)-(30) combined with Eqs. (31)-(32) give a calculable expression 
for the thermal GJT effective potential for quark mass function m{p^). Obviously, it is a function of temperature 
T, quark chemical potential /i, the quark-antiquark condensates (j) and, of course, the current quark mass m as well. 
The result provides a useful means to research chiral symmetry phase transition in QGD in the real-time formalism 
of thermal field theory. 
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